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BL Lac Identification for the Ultraluminous X-ray Source 
Observed in the Direction of NGC 4698 
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Abstract. We report the identification of the optical and radio counterparts of the ultraluminous X-ray (ULX) 
source XMMU J124825. 9+083020 (NGC4698-ULX1). The optical spectrum taken with the VLT yields a redshift 
oi z = 0.43, which implies that the ULX is not associated with the nearby galaxy NGC 4698. The spectral energy 
distribution calculated from the available data indicates that the source is likely to be a BL Lac object. The 
possible synchrotron peak at X-ray energies suggests that this source may be a 7-ray emitter. 
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1. Introduction 

In recent years, the improved imaging capabilities and in- 
creased sensitivity of ROSAT, Chandra and XMM-Newton 
have allowed us to effectively study discrete sources in 
nearby galaxies beyond the Local Group. Particularly in- 
triguing is the discovery of off-nuclear X-ray sources with 
luminosities well above the Eddington limit for a typi- 
cal neutron star, ~ 10'^^ erg s^^, and up to 2 x 10^" erg 
s-i (e.g., Read et al. 1997; Colbert & Mushotzky 1999; 
Roberts & Warwick 2000; Makishima et al. 2000; Fabbiano 
et al. 2001). These sources are typically called ultralumi- 
nous X-ray sources (ULXs). Despite much effort, little is 
presently known about these sources. The identification 
of their optical counterparts is often problematical when 
the sources are superposed against regions of high surface 
brightness in the host galaxy. To date, only two ULXs ap- 
pear to have a clear optical identification (Roberts et al. 
2001; Wu et al. 2002), while for others it has been possible 
to study only the nearby environment (PakuU & Mirioni 
2002; Wang 2002). 

We have started a search for ULXs in a sample of 
nearby galaxies (Foschini et al. 2002). This paper con- 
cerns follow-up observations of the ULX in NGC 4698. We 
discuss the radio and optical counterparts and give a red- 
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shift estimate. We identify the source with a background 
source, most likely a BL Lac object at z = 0.43. 



2. X-ray data (XMM-IMewton) 

NGC 4698 is an Sab spiral galaxy located in the Virgo 
cluster {d = 16.8 Mpc). It hosts an active nucleus, clas- 
sified as a Seyfert 1.9 by Ho et al. (1997). The galaxy 
was observed on 16 December 2001 using the European 
Photon Imaging Camera (EPIC) on board the XMM- 
Newton satellite. EPIC is composed of two instruments: 
the PN-CCD camera (Striider et al. 2001) and two MOS- 
CCD detectors (Turner et al. 2001). The effective exposure 
time was 9.2 ks. 

We find only one ULX apparently associated with 
NGC 4698 (Foschini et al. 2002), located 73" from the 
optical nucleus of the galaxy, at a = 12''48™25.9* and 
5 = -h08°30'20" (J2000). The point source centroid is 
measured to precision of better than O'.'l, and the absolute 
pointing uncertainty is < 4" (Jansen et al. 2001). We de- 
tected 198 photons with MOSl, in a circle of radius 30", 
249 with M0S2, and 414 with PN. These statistics are 
sufficient to perform spectral fitting with simple models. 
The ones we use are the power law (PL), black body (BB), 
bremsstrahlung (BR), and the multicolor black body disk 
(MCD, diskbb in Xspec) by Mitsuda et al. (1984). Fluxes 
were corrected according to the energy encircled fraction 
(Ghizzardi 2001). 
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Fig. 1. i?-band image from VLT-U3/F0RS1, superimposed in contours the smoothed image from the XMM-Newton 
EPIC-M0S2 data in the 0.5 - 10 keV energy band (white over NGC 4698, black in the remaining field). North is up 
and East to the left. The D25 ellipse (3.'8) is shown for comparison. 
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Fig. 2. Archival Hubble Space Telescope image of the re- 
gion around XMMU J124825.9 + 083020 taken in the 
F814W filter. North is up and East to the left. 



For the processing, screening, and analysis of the data 
we used the standard tools in the XMM-SAS software (v. 
5.2) and HEAsoft Xspec (v. 11.0.1). Correction for vi- 
gnetting has not been applied because the source is close 
to the center of the field of view (< 2') and most of the 
detected photons have energies less than 5 keV (see Lumb 
2002). 



Table 1. Results from the fit of the X-ray data. 
Columns: (1) Model: power law (PL), black body (BB), 
bremsstrahlung (BR), multicolor black body disk (MCD); 
(2) free parameter of the model: photon index F for PL, 
temperature (keV) for BB and BR, temperature (keV) 
at the inner disk for MCD; (3) and degrees of free- 
dom of the spectral fitting; (4) flux in the 0.5 — 10 keV 
band (10^^** erg cm~^ s~^); (5) X-ray luminosity in the 
0.5 - 10 keV band (lO''^ erg s'^) calculated for z = 0.43, 
Ho = 75 km s^^ Mpc^^ and qo = 0.5. No absorp- 
tion is required in excess of the the Galactic column, 
Nil = 1-87 X 10^° cm~^. The uncertainties in the param- 
eters are at the 90% confidence level. 



Model 


Parameter 


xVd.o.f. 


Fx 


Lx 


(1) 


(2) 


(3) 


(4) 


(5) 


PL 


F = 2.0 ± 0.2 


18.3/24 


8.6 


3.9 


BB 


fcT = 0.35 ±0.05 


34.5/24 


4.5 




BR 


fcT = 2 ± 1 


19.0/24 


6.7 




MCD 


m„ = 0.6 ±0.1 


23.7/24 


5.4 





The best fit (Fig.|) is found using the power-law model 
with L = 2.0 ± 0.2 (x^ = 18.3, = 24), giving a flux 
of 8.6 X 10~^^ erg cm~^ s~^. If the source is located in 
NGC 4698, for which we assume a distance of 16.8 Mpc, 
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the corresponding luminosity, assuming isotropic emis- 
sion, is ^ 3x 10"^^ erg s^^. For z — 0.43, as discussed below, 
the luminosity becomes 3.9 x 10^^ erg s"^. Additional in- 
formation on the fits using the other models is given in 
Table |ll 




1 2 

channel energy (keV) 



Fig. 3. XMM-Newton spectrum of XMMU J124825.9 + 
083020. Data are from MOSl, M0S2, and PN. 



3. Radio data (VLA) 

The VLA observations of NGC 4698 were performed by 
Ho & Ulvestad (2001) on 29 August 1999 (20 cm) and 
31 October 1999 (6 cm). The source detection at 6 cm 
was reported by Ho & Ulvestad (2001); the flux den- 
sity is 1.13 mJy and the background noise is 0.072 mJy 
beam""*^. The source, which appears largely unresolved at 
a resolution of - 1", is located at a = 12''48'"25.9'' and 
S = -|-08°30'21" (J2000); the radio position is accurate to 
- O'.'l. 

A reanalysis of the VLA data detected the source also 
at 20 cm. It is unresolved with the 1" beam, and its 20 cm 
position is consistent with that measured at 6 cm. We 
find a flux density of 0.34 mJy, with a background noise 
of 0.032 mJy beam^^. The spectral index, defined as Si, oc 
v", is ar « 1.0. 

4. Optical data 

4.1. Digitized Sky Survey (DSS) 

The optical counterpart of the source is visible on the DSS 
and the HST observations of NGC 4698. We are certain 
of the identification from the relative offset between the 
nucleus and the optical counterpart. From the Automatic 
Plate Measurement (APMQ) catalog, the point source has 
R = 18.9 mag and B = 21.9 mag, indicating a very 
red color {B ~ R — 3.0 mag). The original observation 
was made with the 48 inch Schmidt telescope at Palomar 
Observatory on 17 February 1950. 

^ http : //www. ast . cam. ac .uk/'^apmcat/. 



4.2. Hubble Space Telescope (HST) 

We retrieved the cahbrated F450W, F606W, and F814W 
WFPC2 images of NGC 4698 from the HST Archive^ 
to perform photometry of the optical counterpart of the 
source. The observations with the F450W and F814W fil- 
ters were performed on 10 August 2001, while the observa- 
tion with F606W was taken on 31 March 2000. In all the 
filters, the source is located on the WF detector, which 
has a scale of O'.'l per pixel. The total exposure time is 
460 s (2 X 230 s) in the F450W and F814W fihers and 
600 s (400 s and 200 s) in the F606W fiher. We com- 
bined the separate images taken in each filter using the 
crrej task in IRAF to reject cosmic rays. Aperture pho- 
tometry was done using the tasks available in the apphot 
package. The integrated magnitudes, determined using the 
aperture growth curve method, were transformed to the 
Vega magnitude system using the zeropoints given in the 
HST Data. Handbook. We find mF45ow = 21.1 ± 0.1 mag, 
TOpeoew — 19.9 ±0.1 mag, and tof814W — 18.6 ±0.1 mag. 
In the Johnson-Kron-Cousins system these values trans- 
late into B — 21.3 ± 0.1 mag, R = 19.6 ±0.1 mag, and 
/ = 18.5 ±0.1 mag. 

4.3. Las Campanas (Magellan) 

XMMU J124825.9 + 083020 was observed for 900 s 
(3 X 300 s) on 6 April 2002 with the LDSS-2 spectro- 
graph (AUington-Smith et al. 1994) on the 6.5 m Baade 
telescope at Las Campanas Observatory. We observed 
through a 1'.'25 longslit in 0'.'9 seeing, with the slit oriented 
at the parallactic angle. We used the medium red grism, 
which gives a dispersion of 5.3 A pixer^ The FWHM 
spectral resolution was 16 A, and the spatial resolution 
was 0'.'38 pixel^^. 

Basic data reduction was performed using the IRAF 
package. The individual spectroscopic frames were cor- 
rected for overscan, flat-fielded using domeflats and 
summed to obtain the final object frame. The spectral ex- 
traction was done by summing the counts within an aper- 
ture of 6 pixels (2'.'3). Wavelength calibration was achieved 
using a polynomial fit to the lines in the sky spectrum. 

The spectrum in the range 5000 — 7000 A is approxi- 
mately described by the form f\ ~ X^. 

4.4. Very Large Telescope (VLT) 

VLT observations^ were performed with VLT-UT3 
(Melipal) plus the FORSl spectrograph. FORSl is 
equipped with a 2048 x 2048 pixel Tektronix CCD, which 



^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of 
Universities for Research in Astronomy, Inc. under the NASA 
contract NAS 5-26555. 

^ Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 269. D — 5014. 
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Fig. 4. Normalized VLT spectrum with identifications of absorption lines marked. Telluric absorption lines are marked 
with the symbol 0. Shown for comparison is the spectrum of the BL Lac object 0548—322 (Barth et al. 2002), redshifted 
to z = 0.43. 



covers a 6'.8 x 6.'8 field in the standard resolution imaging 
mode with a scale of 0'.'2 pixel"^. 

Two spectra, both with an exposure time of 570 s, were 
acquired on 21 April 2002, starting at 4:52:50 UT. The 
spectra were acquired using Grism #1501 plus order sepa- 
rator GG435, which avoids overlapping of spectral orders 
over a given wavelength; this limited the spectral range 
to 4500 - 9000 A. The slit width was 1" for both spectra, 
and this setup secured a final dispersion of 5.5 A pixel" ^, 
corresponding to a FWHM resolution of 13 A. Before the 
spectroscopic observation, a one minute _R-band acquisi- 
tion image was obtained on the same night starting at 
04:42:52 UT under very good seeing conditions (~0'.'6). 
The object appeared clearly elliptical and extended, with 
a bright core and a fuzzy halo. Moreover, an irregular spot 
located around 1" north of the object core is apparent in 
the VLT image. 

After correction for flat-field and bias, the spectra were 
background subtracted and optimally extracted (Horne 
1986) using IRAF. He-Ne-Ar and Hg-Cd lamps were used 



for wavelength calibration. The wavelength calibration 
was checked against the position of night sky lines; the 
typical error was 0.5 A. Finally, the two spectra were 
stacked together in order to increase the signal-to-noise 
ratio. We encountered problems during observation of the 
spectrophotometric standard, and hence our spectra are 
not flux calibrated. 



Fig. ^ shows the normalized VLT spectrum, with sev- 
eral absorption lines identified [see Laurent-Muehleisen et 
al. (1998) for a discussion on optical identification of BL 
Lacs]. The detected lines include Ca II H&K AA3933, 3968, 
the G band A4304, H/3 A4861, and Fe I A5270: they aU in- 
dicate z « 0.43. The line strengths are heavily diluted by 
the featureless continuum, but the spectrum shares close 
similarity to that of the BL Lac object 0548 — 322 (Barth 
et al. 2002). We shifted the spectrum of 0548 - 322, which 
has a redshift of 0.069, to z = 0.43. 
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Fig. 5. The spectral energy distribution of 
XMMU J124825.9+083020, assembled from data taken 
with the VLA, DSS, HST, and XMM. 



The spectral indices can be used to deduce some gen- 
eral properties of the dominant radiation mechanism. If 
c^x < cti-, the source may exhibit relativistic beaming, 
while if ax > cti- , it may fulfill the conditions of the homo- 
geneous synchrotron model (Harris & Krawczynski 2002) . 
In our case, we have ax ~ ctr ~ 1, so that we cannot 
clearly discriminate between these two cases. If the source 
is a high-frequency peaked BL Lac, however, it is likely 
that the homogeneous synchrotron model is more appli- 
cable. This is confirmed by the axox test of Sambruna et 
al. (1996): in this case, the difference by Oox — Q^x = Q^xox is 
approximately equal to zero, so avoiding a clear discrimi- 
nation between the physical mechanism of the source. 

The spectral energy distribution (Fig. ||) is comparable 
to those of BL Lac objects peaked in the 7-ray domain 
(e.g., Fossati et al. 1998), but the third EGRET catalog 
(Hartman et al. 1999) does not have any source within 
several degrees of XMMU J124825.9 + 083020. 



5. Discussion 

5.1. Classification 

The absorption features detected in the VLT spectrum 
clearly places XMMU J124825.9-h083020 at z = 0.43. We 
identify the source as a BL Lac object, for the following 
reasons. (1) We detect no emission lines; the upper limit 
to any emission feature is ~5 A. (2) The stellar features 
superposed on the featureless continuum have strengths 
generally consistent with those of other BL Lac objects. 
We demonstrate this concretely by comparing a redshifted 
spectrum of the BL Lac object 0548 - 322 (Fig. |). The 
break contrast at 4000 A is ^^29%, again similar to other 
BL Lac objects (Laurent-Muehleisen et al. 1998). (3) The 
optical continuum slope measured in the Magellan spec- 
trum is consistent with those of high-frequency peaked BL 
Lac objects (e.g., Stickel et al. 1993). And (4) the multi- 
wavelength spectral energy distribution (Fig. ||), although 
not assembled from simultaneous observations, is highly 
reminiscent of those of BL Lac objects. 

By using two-point spectral indices, namely the radio- 
to-optical aro and optical-to- X-ray aox, it is possible to 
show that different objects populate different regions of 
the aro — ctox plane (e.g., Brinkmann et al. 1997; Laurent- 
Muehleisen et al. 1999). Another method has been sug- 
gested by Maccacaro et al. (1988), who proposed a nomo- 
graph to link the X-ray flux in the energy band 0.3 — 
3.5 keV and the visual magnitude. The values of aro and 
aox for the present source are 0.42 and 0.95, respectively, 
thus placing it in the region of X-ray selected BL Lacs 
(Brinkmann et al. 1997), or high-energy peaked BL Lacs in 
the arc — ttox diagram of Laurent-Muehleisen et al. (1999). 
The nomograph of Maccacaro et al. (1988) gives a ra- 
tio fx/ fv between 1.3 and 3.8 (depending on whether we 
use the HST optical magnitudes through the F606W or 
F450W filter, respectively), in the regime of AGNs and 
BL Lacs. 



5.2. Search for ULXs and Contamination with 
Background Objects 

This research is a useful demonstration of just how diffi- 
cult it is to identify the physical nature of ULXs. It shows 
the vital importance of redshift determinations. Although 
we had a lot of photometric data on XMMU J124825.9-I- 
083020, by themselves they were insufficient to clearly es- 
tablish whether the source belongs to NGC 4698 or is a 
background object. In the first case, the source would have 
been something similar to a microquasar (e.g., Mirabel & 
Rodriguez 1999), probably located in a globular cluster. 
Indeed, the optical image of the source appeared to be 
slightly extended (angular extent ~ 3"), and its position 
with respect to NGC 4698 suggested that it could be a 
globular cluster, albeit an unusually large one. The pos- 
sibility of an accreting black hole in a globular cluster is 
not so remote: Chandra observations of the globular clus- 
ter system of NGC 4472 show that about 40% of the bright 
low mass X-ray binaries are associated with optically iden- 
tified globular clusters (Kundu et al. 2002). In addition, 
the X-ray luminosity function shows a break near 3 x 10^^ 
erg s~^, suggesting that the brightest X-ray binaries are 
accreting black holes (Kundu et al. 2002), perhaps micro- 
quasars. 

XMM-Newton observations of the Lockman Hole 
(Hasinger et al. 2001) show that the number of background 
sources in the energy band 0.5 — 2 keV with flux greater 
than 4.0 x 10"^** erg cm~^ s~^ (the best fit value from 
Tab. |l]) is 15 deg~^. In the energy band 2 — 10 keV, there 
are 40 sources deg~^ with flux higher than 4.6 x 10""'^"' erg 
cm~^ s~^. Assuming the same log — logS* relation, and 
considering that the D25 area of NGC 4698 is about 
7.9 arcmin^, we expect 0.08 background objects in the 
2 — 10 keV energy band and 0.03 in the energy band 
0.5 — 2 keV. However, despite these low values, we have 
found, in the present case, that the only ULX is a back- 
ground AGN. 
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The above calculations could be underestimated in the 
present case because the VLT images show an unknown 
concentration of galaxies north-east of NGC 4698 (see 
Fig. U), thus suggesting the possibility of a statistically 
meaningful excess of background sources. However, no ad- 
ditional X-ray sources is seen in the present XMM-Newton 
observation. Perhaps, a longer exposure may reveal soft 
X-ray emission or additional sources, if NGC 4698 lies 
along the line of sight to a galaxy cluster. The three X-ray 
sources identified to date, however, have three redshifts: 
XMMU J124825.9 -I- 083020 has z = 0.43, NGC 4698 has 
z = 0.0033, and the 05^ T source IRXS J124828.1 + 
083103 has been recently identified with a Seyfert nucleus 
at z = 0.12 (Xu et al. 2001). Therefore, these three sources 
are not members of a single cluster. 

6. Final Remarks 

We have presented a multiwavelength analysis of the ul- 
traluminous X-ray source XMMU J124825.9-f 083020 ob- 
served in the direction of NGC 4698. The availability of X- 
ray, optical, and radio data makes it possible to construct 
the spectral energy distribution, but the crucial informa- 
tion to decipher the true nature of the source was provided 
by the optical spectrum. We show that the absorption 
features are unambiguously associated with a background 
source at z = 0.43. We argue that the source is most likely 
a BL Lac object, most likely of the high-frequency peaked 
category, based on the absence of emission lines, the di- 
lution of the absorption features, the slope of the optical 
continuum, and the overall shape of the spectral energy 
distribution. 
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